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Abstract 

Q In this paper, achievable rates and optimal resource allocation strategies for imperfectly-known fading relay 
channels are studied. It is assumed that communication starts with the network training phase in which the receivers 
O . estimate the fading coefficients of their respective channels. In the data transmission phase, amplify-and-forward 

and decode-and-forward relaying schemes with different degrees of cooperation are considered, and the corresponding 
achievable rate expressions are obtained. Three resource allocation problems are addressed: 1) power allocation between 

<n : 

Q\ , data and training symbols; 2) time/bandwidth allocation to the relay; 3) power allocation between the source and relay 



o 



X 



in the presence of total power constraints. The achievable rate expressions are employed to identify the optimal 



£C) ■ resource allocation strategies. Finally, energy efficiency is investigated by studying the bit energy requirements in the 

O ' 

qq ' low-SNR regime. 



Index Terms: Relay channel, cooperative transmission, channel estimation, imperfectly-known fading channels, 
achievable rates, optimal resource allocation, energy efficiency in the low-power regime. 

I. Introduction 

In wireless communications, deterioration in performance is experienced due to various impediments 
such as interference, fluctuations in power due to reflections and attenuation, and randomly-varying channel 
conditions caused by mobility and changing environment. Recently, cooperative wireless communications 
has attracted much interest as a technique that can mitigate these degradations and provide higher rates or 
improve the reliability through diversity gains. The relay channel was first introduced by van der Meulen 
in [1], and initial research was primarily conducted to understand the rates achieved in relay channels [2], 
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[3]. More recently, diversity gains of cooperative transmission techniques have been studied in [4]-[7]. 
In [6], several cooperative protocols have been proposed, with amplify-and-forward (AF) and decode-and- 
forward (DF) being the two basic relaying schemes. The performance of these protocols are characterized in 
terms of outage events and outage probabilities. In [8], three different time-division AF and DF cooperative 
protocols with different the degrees of broadcasting and receive collision are studied. In general, the area 
has seen an explosive growth in the number of studies (see e.g., [9], [10], [11], [12], [13], and references 
therein). An excellent review of cooperative strategies from both rate and diversity improvement perspectives 
is provided in [14] in which the impacts of cooperative schemes on device architecture and higher-layer 
wireless networking protocols are also addressed. Recently, a special issue has been dedicated to models, 
theory, and codes for relaying and cooperation in communication networks in [15]. 

As noted above, studies on relaying and cooperation are numerous. However, most work has assumed 
that the channel conditions are perfectly known at the receiver and/or transmitter sides. Especially in mobile 
applications, this assumption is unwarranted as randomly-varying channel conditions can be learned by the 
receivers only imperfectly. Moreover, the performance analysis of cooperative schemes in such scenarios 
is especially interesting and called for because relaying introduces additional channels and hence increases 
uncertainty in the model if the channels are known only imperfectly. Recently, Wang et al. in [16] considered 
pilot-assisted transmission over wireless sensory relay networks, and analyzed scaling laws achieved by the 
amplify-and-forward scheme in the asymptotic regimes of large nodes, large block length, and small snr 
values. In this study, the channel conditions are being learned only by the relay nodes. In [17], Avestimehr 
and Tse studied the outage capacity of slow fading relay channels. They showed that Bursty Amplify-Forward 
strategy achieves the outage capacity in the low SNR and low outage probability regime. Interestingly, they 
further proved that the optimality of Bursty AF is preserved even if the receivers do not have prior knowledge 
of the channels. 

In this paper, we study the achievable rates of imperfectly-known fading relay channels. We assume that 
transmission takes place in two phases: network training phase and data transmission phase. In the network 
training phase, a priori unknown fading coefficients are estimated at the receivers with the assistance of 
pilot symbols. Following the training phase, AF and DF relaying techniques with different degrees of 
cooperation are employed in the data transmission. We first obtain achievable rate expressions for different 
relaying protocols and subsequently identify optimal resource allocation strategies that maximize the rates. 
We consider three types of resource allocation problems: 1) power allocation between data and training 
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symbols; 2) time/bandwidth allocation to the relay; 3) power allocation between the source and relay if 
there is a total power constraint in the system. Finally, we investigate the energy efficiency by finding the 
bit energy requirements in the Iow-snr regime. 

The organization of the rest of the paper is as follows. In Section HH we describe the channel model. 
Network training and data transmission phases are explained in Section [nil We obtain the achievable rate 
expressions in Section [IV] and study the optimal resource allocation strategies in Section [V] We discuss the 
energy efficiency in the Iow-snr regime in Section [VI] Finally, we provide conclusions in Section IVIII 

II. Channel Model 

We consider the three-node relay network which consists of a source, destination, and a relay node. 

Relay 




Source Destination 



Source-destination, source-relay, and relay-destination channels are modeled as Rayleigh block-fading chan- 
nels with fading coefficients denoted by h sr , h sd , and h rd , respectively for each channel. Due to the block- 
fading assumption, the fading coefficients h sr ~ CAf(0, <7 sr 2 ), h sd ~ C7V(0, <7 sd 2 ), and h rd ~ C7V(0, <J rd 2 ^L 
stay constant for a block of m symbols before they assume independent realizations for the following 
block. In this system, the source node tries to send information to the destination node with the help of 
the intermediate relay node. It is assumed that the source, relay, and destination nodes do not have prior 
knowledge of the realizations of the fading coefficients. The transmission is conducted in two phases: network 
training phase in which the fading coefficients are estimated at the receivers, and data transmission phase. 
Overall, the source and relay are subject to the following power constraints in one block: 

||x s , t f + £{||x s || 2 } <mP s , (1) 

||x rit || 2 + £{||x r || 2 }<mP r . (2) 

where x s t and x r t are the source and relay training signal vectors, respectively, and x s and x r are the 
corresponding source and relay data vectors. 

2 x ~ CN{d,a 2 ) is used to denote a proper complex Gaussian random variable with mean d and variance a 2 . 
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III. Network Training and Data Transmission 

A. Network Training Phase 

Each block transmission starts with the training phase. In the first symbol period, source transmits a pilot 
symbol to enable the relay and destination to estimate the channel coefficients h sr and h sd , respectively. In 
the average power limited case, sending a single pilot is optimal because instead of increasing the number of 
pilot symbols, a single pilot with higher power can be used. The signals received by the relay and destination 
are 

y r>t = h sr x Stt + n r , and y d>t = h sd x s , t + n d , (3) 

respectively. Similarly, in the second symbol period, relay transmits a pilot symbol to enable the destination 
to estimate the channel coefficient h rd . The signal received by the destination is 

Vd,t = hrdX r ,t + n r d . (4) 

In the above formulations, n r ~ CAf(0, N ), n d ~ CJ\f(0,N ), and n r d ~ CJ\f(0,N ) represent independent 
Gaussian noise samples at the relay and the destination nodes. 

In the training process, it is assumed that the receivers employ minimum mean- square-error (MMSE) 
estimation. We assume that the source allocates 5 S of its total power for training while the relay allocates 
S r of its total power for training. As described in [25], the MMSE estimate of h sr is given by 

- a 2 rV / 5 s mP s 
a 2 sr d s mP s + N 

where y rjt ~ CJ\f(0 : a 2 r 5 s mP s + N ). We denote by h sr the estimate error which is a zero-mean complex 
Gaussian random variable with variance var(h sr ) = ai . Similarly, for the fading coefficients h sd 

and h rd , we have 

h«= ff / T% -y^ y d ,t~CX(0,a 2 sd 5 s mP s + N ), var(h sd ) = 2 ^ , (6) 
a z sd 5 s mP s + N a z d 6 s mP s + N 

hrd = jl d ^T r M yT ^ ylt-CAf(0,a 2 rd 6 r mP r + N ), var(h rd ) = 2 ^f . (7) 
a^ d d r mP r + N a^ d d r mP r + N 

With these estimates, the fading coefficients can now be expressed as 

h sr hsr ~t~ hsri hsd ^sd ~t~ ^sdi ^rd ^rd ~t~ ^rd- (8) 
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B. Data Transmission Phase 

The practical relay node usually cannot transmit and receive data simultaneously. Thus, we assume that 
the relay works under half-duplex constraint. Hence, the relay first listens and then transmits. As discussed in 
the previous section, within a block of m symbols, the first two symbols are allocated to network training. In 
the remaining duration of m — 2 symbols, data transmission takes place. We introduce the relay transmission 
parameter a and assume that a(m — 2) symbols are allocated for relay transmission. Hence, a can be seen as 
the fraction of total time or bandwidth allocated to the relay. Note that the parameter a enables us to control 
the degree of cooperation. We consider several transmission protocols which can be classified into two 
categories by whether or not the source and relay simultaneously transmits information: non-overlapped and 
overlapped transmission. Note that in both cases, the relay transmits over a duration of a(m — 2) symbols. 
In non-overlapped transmission, source transmits over a duration of (1 — a) (m — 2) symbols and becomes 
silent as the relay transmits. On the other hand, in overlapped transmission, source transmits all the time 
and sends m — 2 symbols in each block. 

We assume that the data vectors x s and x r are composed of independent random variables with equal 
energy. Hence, the covariance matrices of x s are given by 

EU xU — P' I— (1 ~^ )mP * I and EU xU — P' I — (1 ~ 5s)mPs I (9) 

in non-overlapped and overlapped transmissions, respectively. The covariance matrix for x r is 

EU xU = P> I = (1 - 6r)mPr I (10) 
1 r A r (m-2)a 1 } 

1 ) Non-overlapped transmission: We first consider the two simplest cooperative protocols: non-overlapped 
AF, and non-overlapped DF with repetition coding where the relay decodes the message and re-encodes 
it using the same codebook as the source. In these protocols, since the relay either amplifies the received 
signal, or decodes it but uses the same codebook as the source when forwarding, source and relay should 
be allocated equal time slots in the cooperation phase. Therefore, we initially have direct transmission from 
the source to the destination without any aid from the relay over a duration of (1 — 2a) (m — 2) symbols. 
In this phase, source sends x si and the received signal at the destination is given by 

y d i = h sd x s i + n dl . (11) 
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Subsequently, cooperative transmission starts. At first, the source transmits an a (m — 2) -dimensional symbol 
vector x s2 which is received at the the relay and the destination, respectively, as 

y r = /i sr x s2 + n r , and y d2 = h sd x s2 + n d2 . (12) 

For compact representation, we denote the signal received at the destination directly from the source by 
= [yJi yd2\ T where T denotes the transpose operation. Next, the source becomes silent, and the relay 
transmits an a (m — 2) -dimensional symbol vector x r which is generated from the previously received y r [6] 
[7]. This approach corresponds to protocol 2 in [8], which realizes the maximum degrees of broadcasting 
and exhibits no receive collision. The destination receives 

= h rd x r + n r d . (13) 

After substituting the expressions in © into (fTTT)-(fl"3l). we have 

y dl = h sd x sl + h sd x sl + n d i, y r = h sr x s2 + h sr yi s2 + n r , y d2 = h sd x s2 + h sd x s2 + n d2 , (14) 

y r d = h rd K r + h rd x r + n r d . (15) 

We define the source data vector as x s = [x^ xf 2 ] T . Note that we have < a < 1/2 for AF and repetition 
coding DF. Therefore, a = 1/2 models full cooperation while we have noncooperative communications as 
a — > 0. It should also be noted that a should in general be chosen such that a(m — 2) is an integer. 

For non-overlapped transmission, we also consider DF with parallel channel coding, in which the relay 
uses a different codebook to encode the message. In this case, the source and relay do not have to be allocated 
the same duration in the cooperation phase. Therefore, source transmits over a duration of (1 — a)(m — 2) 
symbols while the relay transmits in the remaining duration of a(m — 2) symbols. Clearly, the range of a is 
now < a < 1. In this case, the input-output relations are given by (fT2l) and ([T3l . Since there is no separate 
direct transmission, x s2 = x s and y^ 2 = y^ in (fT2)) . Moreover, the dimensions of the vectors x s , y^, y r are 
now (1 — a)(m — 2), while x r and y^ are vectors of dimension a(m — 2). 

2) Overlapped transmission: In this category, we consider a more general and complicated scenario in 
which the source transmits all the time. In AF and repetition DF, similarly as in the non-overlapped model, 
cooperative transmission takes place in the duration of 2a (m — 2) symbols. The remaining duration of 
(1 — 2a) (m — 2) symbols is allocated to unaided direct transmission from the source to the destination. 
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Again, we have 0<a<l/2in this setting. In these protocols, the input-output relations are expressed as 
follows: 

y dl = h sd x sl + n dl , y T . = h sr x s21 + n r , y d2 = h sd x s21 + n d2 , and y^ = h sd x s22 + h rd x r + n r d . (16) 

Above, x s i, x S 2i, x s22 , which have respective dimensions (1 — 2a) (m— 2), a{m — 2) and a(m — 2), represent 
the source data vectors sent in direct transmission, cooperative transmission when relay is listening, and 
cooperative transmission when relay is transmitting, respectively. Note again that the source transmits all 
the time. x r is the relay's data vector with dimension a{m — 2). y dl , y d2 , y T d are the corresponding received 
vectors at the destination, and y r is the received vector at the relay. The input vector x s now is defined 
as x s = [x^, xj 21 , xf 22 ] T and we again denote y^ = [y dl y d2 \ T . If we express the fading coefficients as 
h = h + h in (fT6l) . we obtain the following input-output relations: 

y d i = h sd x sl + h sd x sl + n d i, y r = h sr x s21 + h sr x s21 + n r , y d2 = h sd x s21 + h sd x s21 + n d2 , and (17) 

y^ = h sd x s22 + h rd x r + h sd x s22 + h rd x r + n d . (18) 

IV. Achievable Rates 

In this section, we provide achievable rate expressions for AF and DF relaying in both non-overlapped 
and overlapped transmission scenarios described in Section HID Achievable rate expressions are obtained by 
considering the estimate errors as additional sources of Gaussian noise. Since Gaussian noise is the worst 
uncorrelated additive noise for a Gaussian model [20], [21], achievable rates given in this section can be 
regarded as worst-case rates. 

We first consider AF relaying scheme. The capacity of the AF relay channel is the maximum mutual 
information between the transmitted signal x s and received signals y d and y^ given the estimates h sr , h sd , h rd : 

C= sup — I(x s ;y d ,y T d \h sr , h sd , h rd ). (19) 

Note that this formulation presupposes that the destination has the knowledge of h sr . Hence, we assume that 
the value of h sr is forwarded reliably from the relay to the destination over low-rate control links. In general, 
solving the optimization problem in (fl9l ) and obtaining the channel capacity is a difficult task. Therefore, we 
concentrate on finding a lower bound on the capacity. A lower bound is obtained by replacing the product 
of the estimate error and the transmitted signal in the input-output relations with the worst-case noise with 
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the same correlation. In non-overlapped transmission, we consider 



z dl = h sd x sl + n dl , z r = h sr x s2 + n r , z d2 = h sd x s2 + n d2 , and z r d = h rd x r + n^, (20) 
as the new noise vectors whose covariance matrices, respectively, are 

E{z dl z dl } = a 2 z \ = a\ s E{^ sl } + N I, E{z r z\} = a 2 z I = ^{x^} + N I, (21) 

E{z d2 z d2 } = a 2 J = al 3 E{x s2 4 2 } + N I, E{z d z r J} = a 2 z A = ^Efax*} + N I. (22) 
Similarly, in overlapped transmission, we define 

Zdl = h sd K s \ + ridi, z r = /i sr x s2 i + n r , z^ 2 = h sd x. s2 i + n^, z^ = h sd x s22 + /i r( iX r + n^, (23) 
as noise vectors with covariance matrices 

E{z dl z ai } = a" z J = al s E{^ sl } + N I, E{z r z\} = <I = ^{x^x^} + N I, (24) 

E{z d2 z d2 } = a 2 Zd I = alE{x s21 4 21 } + N I, E{z d z r J} = a 2 r I = ^^{x^x^J + ^^ix^^ + ^I. 

(25) 

An achievable rate expression is obtained by solving the following optimization problem which requires 
finding the worst-case noise: 



C^I AF = inf sup — J(x a ; y d , y d \h sr , h sd , h rd ). 

Pz dl (-),Pzr(-),Pz d2 (-)>Pzr{-)p Xs (.) Tfl 

The following results provide Iaf for both non-overlapped and overlapped transmission scenarios. 
Theorem 1: An achievable rate of AF in non-overlapped transmission scheme is given by 



(26) 



Iaf = — E 

m 



1 -2a)(m-2)log 1 + 



P' sl \h sd \ 2 



aim — 2) log 1 + 



P' sl \h sd \ 2 



+ f 



Zd2 



P' \h I 2 P'lA I 2 
- r sl|" , sr| - r r |"-rd| 



z d 

(27) 



where 



xy 



P'si\h s d\ 2 _ P'si\h s d\ 2 



1 + x + y 
5 s {l-5 s )m 2 P 2 atJ{l-a) 



z di 



Zd2 



[1 - 5 s )mP s a 2 sd N /(l - a) + (m - 2){a 2 J s mP s + N )N 



\Wsd\ 



(28) 
(29) 
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P' \h 



5 s (l-i> 2 P> s 4 r /(l-« 



P'r\h r d\ 2 



1 - 5 s )mP s a s 2 r iV /(l - a) + (m- 2)(a 2 sr 5 s mP s + N )N 
(5 r (l-^)m 2 P r Vy« 



We 



W r d\ 



(30) 



(3D 



(1 - 5 r )mP r a* d N /a + (m - 2)(afJ r mP r + N )N 

In the above equations and henceforth, w sr ~ £/V(0, 1), w sd ~ CA/"(0, 1), uvd ~ C7V(0, 1) denote indepen 
dent, standard Gaussian random variables. 

Proof: Note that in non-overlapped AF relaying, 



-f( x «; yd, y r d \h sr , h rd ) — i\x*i\ ydi\h S d) + i(x*2] yd2, y d \h sr , h sd , h ra 



(32) 



where the first mutual expression on the right-hand side of (|32l is for the direct transmission and the second 
is for the cooperative transmission. In the direct transmission, we have 



y dl = h sd x. sl + z dl . 



(33) 



In this setting, it is well-known that the worst-case noise z d i is Gaussian [20] and x s i with independent 
Gaussian components achieves 



inf sup I(x al ;y d i|/i sd ) =E 



l-2a)(m-2)log 1 + 



Psi\h s d\ 2 



a 



(34) 



Therefore, we now concentrate on the cooperative phase. For better illustration, we rewrite the channel 
input-output relationships in (fT4l) and (fT5l) for each symbol: 



y r [i] = h sr x s2 [i] + z r [i], y d2 [i] = h sd x s2 [i] + z d2 \i], 



(35) 



for i = 1 + (1 - 2a)(m - 2), (1 - a)(m - 2), and 



y r d [i] = h rd x r [i] + z r d [i\, 



(36) 



for 2 = (1 — a)(m — 2 ) + 1 , . . . , m — 2. In AF, the signals received and transmitted by the relay have following 
relation: 

x r [i] = /3y r [i — aim — 2)], where j3 < \ j — — - _ r " L '7 n ,V J „ r , — — • (37) 



\h sr \ 2 E[\x s2 \ 2 ]+E[\z r W 
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Now, we can write the channel in the vector form 



/ 



y r d [i + a{m - 2)] 



fd.[i] 









H 




x s [i\+ | 




I h r dPh sr J 




A 



1 

h rd f3 1 



B 



Z r [l\ 

V z r d [i + a{m-2)} J 



(38) 



— 



for z = 1 + (1 — 2a) (m — 2), (1 — a)(m — 2), With the above notation, we can write the input-output 
mutual information as 

(l-a)(m-2) 

^(x^sy^^dl^sr^sd^rd) = ^(^[^ydKI^,^,^) = a(m- 2)I(x s ;y d \h sr ,h sd ,h rd ) (39) 

i=l+(l-2a)(m-2) 



where in (|39l) we removed the dependence on i without loss of generality. Note that y is defined in (1381) . 
Now, we can calculate the worst-case capacity by proving that Gaussian distribution for z r , z d2 , and z r d 
provides the worst case. We employ techniques similar to that in [20]. Any set of particular distributions 
for z r , z d2 , and z r d yields an upper bound on the worst case. Let us choose z r , z d 2, and z d to be zero mean 
complex Gaussian distributed. Then as in [6], 

inf sup I(x s ;f d \h sr ,h sd ,h rd ) < Elogdet (I + {E^x.^AA^BE^B*)- 1 ) (40) 

where the expectation is with respect to the fading estimates. To obtain a lower bound, we compute the 
mutual information for the channel in (1381 ) assuming that x s is a zero-mean complex Gaussian with variance 
E(\x s \ 2 ), but the distributions of noise components z r , z d2 , and z r d are arbitrary. In this case, we have 

I\3'8) yrf> \hsri h sd , h rd j h{x s \h sr , h sd , h rd ) h(x s \~y d , h sr , h sd , h rd ) 

^ log7reE(|a; s | 2 ) - log 7re var (x s \y d , h sr , h sd , h rd ) (41) 

where the inequality is due to the fact that Gaussian distribution provides the largest entropy and hence 

K x s\yd, h sr , hd, h rd ) < log Tie var (x s \ y d , h sr , h sd , h rd ). From [20], we know that 



var(x s \y d , h sr , h sd , h rd ) ^ E (x s - x s )(x s - x s y\h sr , h sd , h 



rd 



(42) 



for any estimate x s given y d , h sr , h sd , and h rd . If we substitute the LMMSE estimate x s = R X yR y 1 y d into 
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411) and (|42|) . we obtain 
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I(x s ;y d \h sr ,h s d,hrd)>Elogdet {l+ {EWxs^AA^BE^z^)- 1 ) . 



(43) 



Since the lower bound (1431) applies for any noise distribution, we can easily see that 



inf sup I(x s ; y d \h sr , h sd , h rd ) > Plogdet (I + {E[\x s \ 2 )AA^){BE[zz^)- x ) 

From d40l) and (|44|). we conclude that 



(44) 



inf sup I{x s ; y d \h sr , h sd , h rd ) = Plogdet (I + {E[\x s \ 2 ]AA^){BE[zz^)- 1 ) 



E 



log U + 



1 , -Pjll^sdl 



+ / 



Pil\hsr\ Pr\h r d 



of 



(45) 
• (46) 



In (|46l) . P s i/ and are the powers of source and relay symbols and are given in © and (flOl) . Moreover, 
olt, ,o 2 ,a 2 r are the variances of the noise components defined in (1201) . Now, combining (1261 ), (1321) , (1341) , 
and (|46l ), we obtain the achievable rate expression in (T27T ). Note that (f29T)— (f3T1) are obtained by using the 
expressions for the channel estimates in ©-© and noise variances in (T2TI) and (T22l) . □ 
Theorem 2: An achievable rate of AF in overlapped transmission scheme is given by 



Iaf = — E 

m 



(l-2a)(m-2)log(l + 



P' S 2\h 2 



o 



w| *+(m - 2)alog ( 1 + + /( J ^ 1 ^ 1 



(7 



~rf:2 



a 2 



P s ' 2 |^sd| 2 -P^l^sdl 2 -P^l^rl 2 -P r '|^ 



rd I 



9 ' 2 ' 9 ' 2 

o 0% o Or 

z d2 z i z r Z, 



(47) 



where q(.) is defined as q(a,b,c,d) = "l*^^ . Moreover 



-P^l^dl 2 P.s2\hsd\ 2 



5 s (l-5 s )m 2 P 2 oi d 



a 



~dl 



o 



'-,12 



P's2\k 



(1 - 5 s )mP s o 2 d N + (m - 2){o 2 J s mP s + iV )iVo 
5 s (l-d s )m 2 PX 4 r 



\Wsd\ 



ol (1 - 5 8 )mP 8 alN + (m - 2)(ay s mP s + iV )iV c 



(48) 
(49) 



{m-2){a 2 d S s mP s +N ){a 2 d S r mP r +N )N + {l-S r )mP r a 2 d N Q (a 2 d d s mP s + N )/a + (1 - S s )mP s a 2 d N {a 2 d S r mP r + N Q ) 



3 Here, we use the property that det(I + AB) = det(I + BA) 
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S r (l - S r )m 2 P r 2 a^ 2 sd SsmP s + N )/a\w rd \< 



(m-2)(a 2 d 6 s mP s + N Q ){a 2 d S r mP r + N Q )N + (1 - d r )mP r af d N (a 2 d 5 s mP s + N )/a + (1 - S s )mP s a 2 d N (a 2 d 6 r mP r + N Q ) 

Proof: Note that the only difference between the overlapped and non-overlapped transmissions is that 
source continues its transmission as the relay transmits. As a result, the power of each source symbol is 
now P' s2 given in Q. Additionally, when both the source and relay are transmitting, the received signal at 
the destination is y r d = h sd K s2 2 + Kd^r + h sd K s22 + h rd -x r + n r d . The input-output mutual information in one 
block is 

l( x s;yd,yd\h sr ,h sd ,h rd ) = I{^ s i;ydi\h S d) + I( x s2i,^s22;yd2,y d \h s r,h sd , h rd ). (52) 



The first term on the right-hand-side of (|52l) corresponds to the mutual information of the direct transmis- 
sion and is the same as that in non- overlapped transmission. Hence, the worst-case rate expression obtained 
in the proof of Theorem [Q is valid for this case as well. In the cooperative phase, the input-output relation 
for each symbol can be written in the following matrix form: 



Vd2[i] \ I h sd \ / x s [i] \ I 10 



y d [i + a(m - 2)] / V h rd (3h sr h sd J \ x„[i + a(m — 2)] / \ h rd /3 1 



Z d 2[H 

V z r d [i + a{m-2)] J 



(53) 



where i = 1 + (1 -2a)(m-2), (1 - a)(m - 2) and (3 ^ </ .- I2 J [ BS— Note that we have defined 

V \hsr\ E[\x s \ \-\-E[\z r \ \ 

x s = [xj 1? xj 21 , xJ 22 ] T , and the expression in (1531) uses the property that x 2 i(j) = x s (j + (1 — 2a) (m — 2)) 
and x S 22(j) = x s (j + (1 — a)(m — 2)) for j = 1, . . . , a(m — 2). The input-output mutual information in the 
cooperative phase can now be expressed as 

(l-a)(m-2) 

Xs2i,Xs22;yd2,y r d,\h S r,h sd ,h rd ) = ^(x s [i];y rf [i]|/i sr ,/i S(i ,/ir<i) = ai(m-2)I(x s ;y d \h sr h sd ,h rd ) (54) 

i=l+(l-2a)(m-2) 



where in (1541) we removed the dependence on i without loss of generality. Note that x and y are defined in 
(1531) . As shown in proof of Theorem [H the worst-case achievable rate for cooperative transmission is 

inf sup I{* s -y d \h sr XdXd) = £logdet (I + {E[± s ±l]AA^){BE[2,z^y l ) . (55) 
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Using the definitions in (1531) and evaluating the logdet expression in (1551) . and combining the direct 
transmission worst-case achievable rate, we arrive to (l47l ). (l48l)-(l5T|) are obtained by using the expressions 
for the channel estimates in ©-([7]) and noise variances in (1241) and (l25l ). □ 
Next, we consider DF relaying scheme. In DF, there are two different coding approaches [7], namely 
repetition coding and parallel channel coding. We first consider repetition channel coding scheme. The 
following results provide achievable rate expressions in both non-overlapped and overlapped transmission 
scenarios. 

Theorem 3: An achievable rate expression for DF with repetition channel coding for non-overlapped 
transmission scheme is given by 



'DFr 



(l-2a)(m-2) 



E 



m 



log 1 + 



P' s i\hsd\ 



a 



-dl 



aim — 2) 
+ — J -mm{h,I 2 } 



m 



where 



h = E 



P' \h 

r sl\ IL s 



log ( 1 + ^ 



and In — E 



log (l 



P' sl \hJ 2 PL\k 12 



+ 



a 



~<V1 



(56) 



(57) 



Moreover, p 'si\^\\ p 's^\\ p iilM 2 ? and K\K£ are the same as defined in ^5J_^TJ. 

Proof: For DF with repetition coding in non-overlapped transmission, an achievable rate expression is 



I( x si, ydi\h s d) + min <^ 7(x s2 ; y r \h sr ), ^(x s2 ; Yd, y r d\h S d, h 



rdi 



(58) 



Note that the first and second mutual information expressions in (1581) are for the direct transmission between 
the source and destination, and direct transmission between the source and relay, respectively. Therefore, as 
in the proof of Theorem [Q the worst-case achievable rates can be immediately seen to be equal to the first 
term on the right-hand side of (|56l) and J l5 respectively. 

In repetition coding, after successfully decoding the source information, the relay transmits the same 
codeword as the source. As a result, the input-output relation in the cooperative phase can be expressed as 



Vd[i\ 

y r d [i + a(m - 2)] 









H 




x,[i] + 




V h rd {3 j 





Zd2[l>\ 

z r d [i + a(m-2)] 



(59) 



ydW 



— V" 

z\i\ 



where (3 < J ^p^pj- From (l59l , it is clear that the knowledge of h sr is not required at the destination. We 
can easily see that (l59l) is a simpler expression than (|38l ) in the AF case, therefore we can adopt the same 
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methods as employed in the proof of Theorem \T\ to show that Gaussian noise is the worst noise and I2 is 
the worst-case rate. □ 
Theorem 4: An achievable rate expression for DF with repetition channel coding for overlapped trans- 
mission scheme is given by 

P's2\hsd\ 2 



(l-2a)(m-2) E1 

1 DFr — & 



m 



log 1 + 



a 



+ J —mm{hJ 2 } 

m 



(60) 



where 



h = E 



log 1 + 



P' \h I 2 

1 s2\"'sr\ 



h = E 



log 1 + 



P's2^ l sd\ . Pr\h r d\ . P^l^sdl . P^l^sd] P&V 1 



+ 



+ 



a 



~d:l 



a,. 



air 



+ 



sd I 



Zd2 Z d 



(61) 



P's2\hsd\ P' s2 \h sd \ ^ P' s2 \h sr \ P' s2 \h sd \ ^ P'\h rd \ 2 haye ^ same expressions as m ggj_gj]). 

CT -5 a *d 

Proof: Note that in overlapped transmission, source transmits over the entire duration of (m — 2) symbols, 
and hence the channel input-output relation in the cooperative phase is expressed as follows: 



VdXH 

y r d [i + a{m - 2)] 



h s d 
hrdft h sd 



X s [l\ 

x 8 [i + a(m — 2)} 



zmm 
z r d2 [i + aim - 2)} 



(62) 



— 

X \i] 



— V" 

z\i] 



The result is immediately obtained using the same techniques as in the proof of Theorem |2l □ 
Finally, we consider DF with parallel channel coding and assume that non-overlapped transmission scheme 
is adopted. From [11], we note that an achievable rate expression is 

min{(l - a)/(x a ; y r \h sr ), (1 - a)/(x a ; y d \h s d) + crf(x r ; y r d \h rd )}. 



Note that we do not have separate direct transmission in this relaying scheme. Using similar methods as 
before, we obtain the following result. The proof is omitted to avoid repetition. 

Theorem 5: An achievable rate of non-overlapped DF with parallel channel coding scheme is given by 



(l-q)(m-2) 

1df p — mm < h, 

' m 



log 1 + 



P' \h 



[l-a)(m-2) 



E 



m 



log 1 + 



P' sl \h sd \ 2 



a 



~d2 



a(m — 2) 



E 



rn 



log 1 



P'r\h rd \ 2 



air 



+ 



(63) 
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where p °^* d] \ p '^\\ an( i ^ g i ve n in (HUl-dlB- □ 

z d2 "r 

V. Optimal Resource Allocation 

Having obtained achievable rate expressions in Section [TV] we now identify optimal resource allocation 
strategies that maximize the rates. We consider three resource allocation problems: 1) power allocation 
between the training and data symbols; 2) time/bandwidth allocation to the relay; 3) power allocation 
between the source and relay under a total power constraint. 

We first study how much power should be allocated for channel training. In nonoverlapped AF, it can 
be seen that 5 r appears only in Pr j^ in the achievable rate expression (l27l) . Since f(x,y) = x ^ +y is 

. . P'\h 2 

a monotonically increasing function of y for fixed x, (1271) is maximized by maximizing rl 2 rdl . We can 
maximize ^^jH by maximizing the coefficient of the random variable \w rc i\ 2 in (|3T| ), and the optimal S r 
is given below: 



&opt _ ~mP r (j 2 rd ~ amNp + 2ajV + ygm - 2)(m 2 P r (J 2 rd aN + m?P?o* d + amNg + mP r a 2 d N - 2mP r a 2 d aN - 2N^a) 

mP r a 2 d ( — 1 + am — 2a) 

Optimizing S s is more complicated as it is related to all the terms in (|27T ), and hence obtaining an analytical 

P 1 \h j I 2 P' \h I 2 

solution is unlikely. A suboptimal solution is to maximize sl 2 and al 2 separately, and obtain two 

z dl azr 

solutions o s l and o s2 , respectively. Note that expressions for o sl and o s2 are exactly the same 
as that in (|64|) with P r , a, and a r d replaced by P s , (1 — a), and a s d and a sr , respectively. When the 
source-relay channel is better than the source-destination channel and the fraction of time over which direct 
transmission is performed is small, Pj,1 ^ s ^ is a more dominant factor and <5™ &opt i s a good choice for training 
power allocation. Otherwise, $ s s u ^ opt might be preferred. Note that in non-overlapped DF with repetition and 

P'\h I 2 

parallel coding, rl 2 ' 1 is the only term that includes o r . Therefore, similar results and discussions apply. For 
instance, the optimal 5 r has the same expression as that in (l64l) . Figure Q] plots the optimal 5 r as a function 
of a rd for different relay power constraints P r when m = 50 and a = 0.5. It is observed in all cases 
that the allocated training power monotonically decreases with improving channel quality and converges to 
^am-2a-i ~ ~ 0.169 which is independent of P r . 

In overlapped transmission schemes, both 5 S and 5 r appear in more than one term in the achievable rate 
expressions. Therefore, we resort to numerical results to identify the optimal values. Figures |2] and [3] plot 
the achievable rates as a function of S s and 5 r for overlapped AF. In both figures, we have assumed that 
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a s d = l,a sr = 2,a r d = 1 and m = 50, iVo = 1, a = 0.5. While Fig. [2l where P s = 50 and P r = 50, 
considers high snrs, we assume that P s = 0.5 and P r = 0.5 in Fig. [3j In Fig. [2l we observe that increasing 
S s will increase achievable rate until 5 S ~ 0.1. Further increase in 5 S decreases the achievable rates. On 
the other hand, rates always increase with increasing S r . This indicates that cooperation is not beneficial in 
terms of achievable rates and direct transmission should be preferred. On the other hand, in the low-power 
regime considered in Fig. [3l the optimal values of 5 S and 5 r are approximately 0.18 and 0.32, respectively. 
Hence, the relay in this case helps to improve the rates. 

Next, we analyze the effect of the degree of cooperation on the performance in AF and repetition DF. 
Figures |4]|7] plot the achievable rates as a function of a which gives the fraction of total time/bandwidth 
allocated to the relay. Achievable rates are obtained for different channel qualities given by the standard 
deviations a sdl a sr , and a r( i of the fading coefficients. We observe that if the input power is high, a should 
be either 0.5 or close to zero depending on the channel qualities. On the other hand, a = 0.5 always gives 
us the best performance at low snr levels regardless of the channel qualities. Hence, while cooperation is 
beneficial in the Iow-snr regime, noncooperative transmissions might be optimal at high snrs. We note from 
Fig. |4] that cooperation starts being useful as the source-relay channel variance a 2 sr increases. Similar results 
are also observed in Fig |5] Hence, the source-relay channel quality is one of the key factors in determining 
the usefulness of cooperation in the high snr regime. 

In Fig. [8l we plot the achievable rates of DF parallel channel coding, derived in Theorem |5] We can see 
from the figure that the best performance is obtained when the source-relay channel quality is high (i.e., 
when <j s d = 1, a sr = 10, a r a = 2). Additionally, we observe that as the source-relay channel improves, 
more resources need to be allocated to the relay to achieve the best performance. We note that significant 
improvements with respect to direct transmission (i.e., the case in which a — > 0) are obtained. Finally, we 
can see that when compared to AF and DF with repetition coding, DF with parallel channel coding achieves 
higher rates. On the other hand, AF and repetition coding DF have advantages in the implementation. 
Obviously, the relay, which amplifies and forwards, has a simpler task than that which decodes and forwards. 
Moreover, as pointed out in [14], if AF or repetition coding DF is employed in the system, the architecture 
of the destination node is simplified because the data arriving from the source and relay can be combined 
rather than stored separately. 

In certain cases, source and relay are subject to a total power constraint. Here, we introduce the power 
allocation coefficient 6, and total power constraint P. P s and P r have the following relations: P s = OP, 
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P r = (1 — 9)P, and P s + P r < P. Next, we investigate how different values of 9, and hence different power 
allocation strategies, affect the achievable rates. An analytical results for 9 that maximizes the achievable 
rates is difficult to obtain. Therefore, we again resort to numerical analysis. In all numerical results, we 
assume that a = 0.5 which provides the maximum of degree of cooperation. First, we consider the AF. 
The fixed parameters we choose are P = 100, N — 1, 5 S — 0.1, 5 r = 0.1. Fig. [9] plots the achievable rates 
in the overlapped transmission scenario as a function of 9 for different channel conditions, i.e., different 
values of a sr ,a r d, and a s d- We observe that the best performance is achieved as 9 — > 1. Hence, even in 
the overlapped scenario, all the power should be allocated to the source and direct transmission should be 
preferred at these high SNR levels. Note that if direct transmission is performed, there is no need to learn the 
relay-destination channel. Since the time allocated to the training for this channel should be allocated to data 
transmission, the real rate of direct transmission is slightly higher than the point that the cooperative rates 
converge as 9 — > 1. For this reason, we also provide the direct transmission rate separately in Fig. [9] Further 
numerical analysis has indicated that direct transmission over performs non-overlapped AF, overlapped and 
non-overlapped DF with repetition coding as well at this level of input power. On the other hand, in Fig. [10] 
which plots the achievable rates of non-overlapped DF with parallel coding as a function of 9, we observe 
that direct transmission rate, which is the same as that given in Fig. |9] is exceeded if a sr = 10 and hence 
the source-relay channel is very strong. The best performance is achieved when 9 « 0.7 and therefore 70% 
of the power is allocated to the source. 

Figs. [FQ H2l and [13] plot the non-overlapped achievable rates when P = 1. In all cases, we observe 
that performance levels higher than that of direct transmission are achieved unless the qualities of the 
source-relay and relay-destination channels are comparable to that of the source-destination channel (e.g., 
a s d = 1, a sr = 2, a rd = 1). Moreover, we note that the best performances are attained when the source- 
relay and relay-destination channels are both considerably better than the source-destination channel (i.e., 
when a s d = 1, a sr = 4, a r d = 4). As expected, highest gains are obtained with parallel coding DF although 
repetition coding incur only small losses. Finally, Fig. [14] plot the achievable rates of overlapped AF when 
P — 1. Similar conclusions apply also here. However, it is interesting to note that overlapped AF rates are 
smaller than those achieved by non-overlapped AF. This behavior is also observed when DF with repetition 
coding is considered. Note that in non-overlapped transmission, source transmits in a shorter duration of time 
with higher power. This signaling scheme provides better performance as expected because it is well-known 
that flash signaling achieves the capacity in the low-SNR regime in imperfectly known channels [18]. 
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VI. Energy efficiency 

Our analysis has shown that cooperative relaying is generally beneficial in the low-power regime, resulting 
in higher achievable rates when compared to direct transmission. In this section, we provide an energy 
efficiency perspective and remark that care should also be taken when operating at very low snr values. The 
least amount of energy required to send one information bit reliably is given byD || = C (*snr) where C (snr) 
is the channel capacity in bits/symbol. In our setting, the capacity will be replaced by the achievable rate 
expressions and hence the resulting bit energy, denoted by -j^-, provides the least amount of normalized 
bit energy values in the worst-case scenario and also serves as an upper bound on the achievable bit energy 
levels in the channel. 

We note that in finding the bit energy values, we assume that snr = P/N where P = P r + P s is the 
total power. The next result provides the asymptotic behavior of the bit energy as snr decreases to zero. 

Theorem 6: The normalized bit energy in all relaying schemes grows without bound as the signal-to-noise 
ratio decreases to zero, i.e., 



E, 



b,U 



V SNR 1 

= lim — = = oo. (65) 

7=0 SNR ^° ASNR) 1(0) 



iV, 

Proof: The key point to prove this theorem is to show that when snr — > 0, the mutual information decreases 
as snr 2 , and hence 1(0) = 0. This can be easily shown because when P — > 0, in all the terms Psl ^ sd ^ > j^ijfeg! ; 

a *&\ °" z d2 

P'Ahsr\ 2 P'\h rd \ 2 P' 2 \h 3d \ 2 PL\h sd \ 2 PU\h sr \ 2 P'Sh sd \ 2 , P'\h rd \ 2 ■ rj,, rnr^] , • . 

— 1_ _ Ll _r d i — , 3 i sa\ — s2 \ sa t _s2i_ — 1_ ? s2l^ ? an( j t rdi m Theorems QJ|5] the denominator goes to a 



T 2 ) ^3 ' „2 J '2 

z dl z d2 



-d 

constant while the numerator decreases as P 2 . Hence, these terms diminish as snr 2 . Since log(l + x) = 
x + o(x) for small x, we conclude that the achievable rate expressions also decrease as snr 2 as snr vanishes. 

□ 

Theorem [6] indicates that it is extremely energy-inefficient to operate at very low snr values. We identify 
the most energy-efficient operating points in numerical results. We choose the following numerical values 
for the fixed parameters: 5 S = 5 r = 0.1, a sd = 1, o sr = 4, oy rf = 4, a = 0.5, and 9 = 0.6. Fig. [i~5 
plots the bit energy curves as a function of snr for different values of m in the non-overlapped AF case. 
We can see from the figure that the minimum bit energy, which is achieved at a nonzero value of snr, 
decreases with increasing m and is achieved at a lower snr value. Fig. [16] shows the minimum bit energy 
for different relaying schemes with overlapped or non-overlapped transmission techniques. We observe that 

4 Note that jj^ is the bit energy normalized by the noise power spectral level No. 



18 



the minimum bit energy decreases with increasing m in all cases . We realize that DF is in general much 
more energy-efficient than AF. Moreover, we note that employing non-overlapped rather than overlapped 
transmission improves the energy efficiency. We further remark that the performances of non-overlapped DF 
with repetition coding and parallel coding are very close. 

VII. Conclusion 

In this paper, we have studied the imperfectly-known fading relay channels. We have assumed that the 
source-destination, source-relay, and relay-destination channels are not known by the corresponding receivers 
a priori, and transmission starts with the training phase in which the channel fading coefficients are learned 
with the assistance of pilot symbols, albeit imperfectly. Hence, in this setting, relaying increases the channel 
uncertainty in the system, and there is increased estimation cost associated with cooperation. We have 
investigated the performance of relaying by obtaining achievable rates for AF and DF relaying schemes. We 
have considered both non-overlapped and overlapped transmission scenarios. We have controlled the degree 
of cooperation by varying the parameter a. We have identified the optimal resource allocation strategies 
using the achievable rate expressions. We have observed that if the source-relay channel quality is low, then 
cooperation is not beneficial and direct transmission should be preferred at high snrs. On the other hand, 
we have seen that relaying generally improves the performance at low snrs. We have noted that DF with 
parallel coding provides the highest rates. Additionally, under total power constraints, we have identified the 
optimal power allocation between the source and relay. We have again pointed out that relaying degrades 
the performance at high snrs unless DF with parallel channel coding is used and the source-relay channel 
quality is high. The benefits of relaying is again demonstrated at low snrs. We have noted that non-overlapped 
transmission is superior compared to overlapped one in this regime. Finally, we have considered the energy 
efficiency in the low-power regime, and proved that the bit energy increases without bound as snr diminishes. 
Hence, operation at very low snr levels should be avoided. From the energy efficiency perspective, we have 
again observed that non-overlapped transmission provides better performance than overlapped transmission. 
We have also noted that DF is more energy efficient than AF. 
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Fig. 1. 5 r vs. a r d for different values of P r when m — 50. 




Fig. 2. Overlapped AF achievable rates vs. 8 S and S r when P s = P r — 50 
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Fig. 3. Overlapped AF achievable rates vs. 5 S and 8 r when P s = P r = 0.5 




Fig. 4. Overlapped AF achievable rate vs. a when P s — P r — 50, 8 S — S r — 0.1, m = 50. 
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Fig. 5. Overlapped DF with repetition coding achievable rate vs. a when P 3 = P r = 50, S s = S r = 0.1, m = 50. 
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Fig. 6. Overlapped AF achievable rate vs. a when P s = P r — 0.5, c5 s = Sr = 0.1, m = 50. 
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Fig. 7. Overlapped DF with repetition coding achievable rate vs. a when P 3 = P r = 0.5, 5 S = S T = 0.1, m = 50. 




Fig. 8. Non-overlapped DF parallel coding achievable rate vs. a when P s = P r = 0.5, S s = 5 r = 0.1, m = 50. 



24 




Fig. 9. Overlapped AF achievable rate vs. 9. P = 100, m = 50. 




Fig. 10. Non-overlapped Parallel coding DF rate vs. 9. P = 100, m = 50. 
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Fig. 11. Non-overlapped AF achievable rate vs. 9. P = 1, m = 50. 




Fig. 12. Non-overlapped Repetition coding DF rate vs. 0. P = 1, m = 50. 
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Fig. 13. Non-overlapped Parallel coding DF rate vs. 9. P = 1, m = 50. 




Fig. 14. Overlapped AF achievable rate vs. 9. P = 1, m = 50. 
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Fig. 15. Non-overlapped AF Eb,u/No vs. SNR 




Fig. 16. Eb,u/No vs. m for different transmission scheme 
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